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Jenkins has pointed out that the process t — > cW + W~ is GIM suppressed in the 
standard model. In this note, we calculate the branching ratio for a wide range of 
models, in which the decay occurs at tree level through exchange of a scalar, fermion or 
vector. In the case of scalar exchange, a scalar mass between 2mw and 200 GeV leads to 
a resonant enhancement, giving a branching ratio as high as a few tenths of a percent. 
We then note that all of these models will also allow W + W~ — > ic + tc, and we calculate 
the single-top/single-charm production rate at the LHC. The rates are such that the 
background from single-top/single-bottom production will probably swamp the signal. 
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1 Introduction 



It has now been established that the mass of the top quark is well above 163 GeV. This 
makes the decay t — > c W + W~ kinematically accessible. The branching ratio for this 
decay in the Standard Model was studied by JenkinsjTJJ. The diagram is given in Figure 
1, where the internal fermion can be any of the charge —1/3 quarks. Jenkins showed that 
the rate is extremely small. There are two reasons for this. The first, of course, is phase 
space. The second is a GIM cancellation, which makes the resulting rate proportional 
to the square of the 6-quark mass. Jenkins noted that this second suppression might be 
absent in other models, and that in these models the rate might be measurable. 

In this letter, we first look at the decay t — > c W + W~ in a very generic set of models, 
in which the exchanged particle can be either a scalar, a fermion or a vector boson. We 
then note that in any model in which the decay t — > c W + W~ occurs, one will also 
have the process W + W~ — > ci + tc. This process will not be phase-space suppressed, 
and could thus be much larger. We will calculate the rate for p p — > W + W~ — *> ci +tc 
for the same set of models.Q 

In section 2, we discuss the generic set of models that will be considered. The rate 
for £ — > c W + W~ in these models will be calculated in section 3, and the rate for 
p p — > W + W~ — > ci+tc will be determined in section 4. Section 5 contains a discussion 
of the results and conclusions. 



2 Models 

The most general set of operators involving a t, c, W + and W~ will have 36 coefficients, 
each an arbitrary function of the appropriate Mandelstam variables. Rather than try 
to deal with this large set, we will restrict the discussion to models in which the t — > 
c W + W~ and W + W~ —>■ ci+tc processes can occur at tree-level. This seems reasonable; 
the rates are already very small, and in models in which these processes arise from loops 
and higher- dimensional operators, one would expect the rates to be even smaller. Thus, 
for example, the coupling of a fermion to the W^-bosons will be taken to be the most 
general combination of V and A but not T, since tensor couplings generally arise from 
loops. Other than that restriction, the models will be completely general; the exchanged 
particle can either be a scalar, fermion or vector boson. 

One of the simplest extensions of the standard model is the two-Higgs model, in which 
an additional scalar doublet is included. In general, such an extension will have tree-level 
flavor-changing neutral currents. One can avoid such currents by either coupling all 
of fermions to one doublet, or by coupling the charge 2/3 quarks to one doublet and the 
charge —1/3 quarks to the other. This can be done via a discrete symmetry. However, it 
has been pointed out|4| that there is no phenomenological need to do this, as long as the 

1 The related process e + e~ — > ci + tc + vv as well as the rate for t — > c W + W~ , in the particular 
case of scalar exchange, has been calculated by Bar-Shalom et al.@] 
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couplings of the scalars to the light quarks is small. From an analysis of mass matrices, 
Cheng and Slier [|J argued that the most natural value for the flavor-changing couplings 
of the neutral scalars to fermions is the geometric mean of the Yukawa couplings of the 
fermions, i.e. the top-charm-scalar coupling should be of 0(J (gY)t(gy)c), where gy is the 
conventional Yukawa coupling. From this ansatz, one finds that the lower bound on the 
exchanged scalar mass is fairly weak, of 0(100 — 300) GeV (the precise bound depends on 
the amount of fine-tuning one is willing to tolerate). We will scale the top-charm-scalar 
coupling constant by this factor. 

The diagram for scalar exchange is given in Figure 2. We have left the coupling 
arbitrary. As discussed above, the most natural value for a and b is 0(1), although in 
principle they could be larger (perturbation theory would become questionable if they 
were much larger than O(10)). Since the H.iggs-W + -W~ coupling depends on the Higgs 
vacuum expectation value, one expects C w to be a ratio of vacuum values, i.e. of 0(1). 

The diagram for fermion exchange is given in Figure 1. In the standard model, we 
expect a = a' = — 1, and the contributions of the d,s and b quarks all cancel in the 
massless limit, leading to a result that the diagram is GIM suppressed. If there is a 
fourth generation, however, the internal fermion must be heavy. 

The diagram for vector exchange is given in Figure 3. Here, we expect the contribution 
to be small. In order to interact with two W's, the neutral vector boson must mix with 
the Z-boson. Since the properties of the Z are in stunning agreement with theoretical 
expectations, any such mixing must be very small. In a wide range of models||, the 
mixing must be much less than 0.01, and thus we expect the effective coupling, g 1 , to be 
less than one percent of the weak gauge coupling. The coupling g 2 can be large, but the 
mass of the vector boson must be greater than 170 GeV to avoid dominating top decays. 

We now turn to the calculation of t — > c W + W~ in these models. 



3 t -> c W + W~ 

The top quark mass is known to be between (approximately) 170 and 180 GeV. In order 
to get the largest plausible value for the branching ratio for t — > c W + W~ , we will take 
it to be 180 GeV. The differential decay rate is given by 

dT \M\ 2 



dE + dE_ 647r 3 M t 

where E± are the energies of the emitted W^. Dividing by the standard model decay 
rate of the top quark, which to within a few percent is given by G F M?/(8ny/2), gives the 
branching ratio, after integrating over E±. For each of the models discussed above, we 
will give the square of the matrix element, and plot the resulting branching ratio. The 
polarizations of the final state Ws will be summed over, and the charmed quark mass 
will be set equal to zero, except when it appears as a parameter in the scalar Yukawa 
coupling. 
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For the model with scalar exchange, shown in Figure 2, the matrix element squared 
is given by 

g 4 M t M c C 2 w (a 2 + b 2 ) 



11 ~ (s - M 2 ) 2 + r|M s 2 P * ' Pc 
and s = M t (2E^ 



2m 2 w 



-1+2 



(2) 



where p t -p c = M t (M t -E+- EJ) and s = M t (2E + + 2E_M t ). For the scalar width, we 
use the expression of Atwood et al.||, and include the effects of S — > WW* (which has 
a small but noticeable effect for M s very near 2mw) Inserting into the expression for the 
branching ratio gives the result in Figure 4. There, we have taken C^ia 2 + b 2 ) = 1. 

One can see the resonance where the scalar mass is between 2m w and M t . In this 
region, the scalar can be on-shell, leading to a much larger rate. In this case, the decay 
would be clearly measureable. Outside of the resonance region, however, the rate falls 
very rapidly.0 

The diagram for fermion exchange is given in Figure 1. The minimum value for the 
fermion mass, Mf, is 130 GeV. If we define a = 1 + aa'; b = a + a'; a' = 1 — ««' and 
b' = a — a', then the matrix element squared can be written as 



\M\ l 



9i9i\Vcf\ 2 \Vft\ 
A(k 2 - Mj) 2 



(M 2 (a 2 + b 2 )T x + M 2 (a' 2 + b' 2 )T 2 - M f M t (aa' + bb')T 3 ) (3) 



where k 2 = M 2 + m 2 ^ — 2M t E + . Note that in the standard model, a 



a 



-1, so 



a = 
are 



-6 = 2 and a' = b' = 0, thus the T 2 and T3 contributions vanish. The contributions 



T — 1 

11 ~ Aw 2 



where w 



m 



4w 2 + 5w — 3w 3 — 5xw 2 + 3x — llxw + 3w 3 x — 16x 2 
—Ax 2 w + 28x 3 + 12x 3 w — lQx 4 — 2yw 2 — 5yw + 3w 3 y 
—2xyw 2 — 2xy + \2xyw + 8x 2 y — Ax 2 yw — 8x 3 y 

x = E + , y = E- and we have chosen units where M t = 1, 



(4) 



1 

Aw 2 



Ax 3 



and 



Aw 2 



w — w x — x — wx + Ax + Awx' 
—w 2 y + wy + 2xy — Awxy — Ax 2 y 



llw 2 — 3w + 26w 2 x + 12wx — 12wx 2 



(5) 



Aw 2 y 



(6) 



The results for the branching ratio, where we have set g\ and g 2 to be the electroweak 
coupling, and \V c f \ = \Vft\ = 1 are given in Figure 5. The solid curve is the contribution 

2 This resonance region was also noted in the work of Bar-Shalom, et al.0. Our result for this decay 
is in agreement with theirs. 
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from T\ (choosing a 2 + b 2 = 1), the dashed curve is the contribution from T2 (choosing 
a' 2 + b' 2 = 1) and the dotted curve is the contribution from T3 (choosing aa! + bb' = 1). 
Note that T 3 is defined so that the contribution to the square of the matrix element is 
negative; it is an interference term and we have verified that the total rate is always 
positive. Although we have assumed that the V c f and Vft KM matrix elements are 
unity (and the cross section obviously scales as the square of each) , they are constrained 
somewhat by the unitarity of the KM matrix. These constraints are quite weak, giving 
\V c f \ < 0.6 and Vft is essentially unconstrained. (Should the value of \V c f \ be a small as 
one generally expects for mixing across two generations, then the branching ratio would 
be several orders of magnitude smaller.) Note that the resulting branching ratio is very 
small, generally much less than a part in 10 5 . For light masses, the branching ratio 
increases as one approaches the resonance region (between m c + my/ and m t — myy ), but 
the current lower bound on the mass eliminates the possibility that the heavy fermion 
is in that region] 3 ] Thus, in the case of fermion exchange, it appears that the process 
t — > c W /+ W /_ will not be measureable, unless \V c f\ is surprisingly large and the heavy 
fermion has a mass very near its current lower bound. 

For vector exchange, the diagram is given in Figure 3. The expression for the square 
of the matrix element is 

l-^l 2 = L^m^V - i - 12^ 3 - lW + 7 f+ 7 -f + 12xw 3 

+2Ax 2 w 2 - Ax 2 + lAx 2 w + 3x 3 - 12x 3 w - I9yw 2 

—Axw — Ayw + 12w 3 y + 2Axyw 2 — lOxy + 20xyw 

+17x 2 y - 20x 2 yw - 8x 3 y + 2Ay 2 w 2 - Ay 2 + lAy 2 w 

+17xy 2 — 20xy 2 w — 16x 2 y 2 + 3y 3 — 12y 3 w — 8xy 3 } 

(7) 

where w = rriyy, x = E + , y = E_ and we have chosen units where M t = 1. Putting this 
expression into the branching ratio gives the result in Figure 6. Here, we have expressed 
the results obtained by setting g 2 g 2 (l + a 2 ) = 1. As discussed in section 2, however, 
one expects g 2 to be two orders of magnitude less than the weak gauge coupling, so as 
not to significantly mix with the Z. Thus, g\ is at most 10~ 5 , and thus the decay is 
unmeasureable . 

Part of the reason that these decay rates are so small is simply phase space. Any 
model with t — > c W + W~ will also have W + W~ —>■ ci+tc and one could look for single- 
top/single-charm production in hadron colliders. Such a process will not be phase-space 
suppressed and might be much larger. We now turn to this process. 

3 D0 has quoted]?! a firm lower bound of mz + m\, on the mass; from their data it appears that a 2a 
lower bound would be approximately 110 — 120 GeV. 
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4 W + W~ ct + tc 

In order to relate a process involving W + W~ scattering to a process involving scatter- 
ing, one needs to know the W content of the proton. We follow the procedure of Johnson, 
Olness and Tung.|| The cross section is 

a{p P - WW) = J d yi dy 2 d Zl dz 2 fU^fAyJfwi^f^Hzi^yms) (8) 

where, for example, fyy is the W content of the down quark. Changing variables to 
x\ = yiZ\, X2 = 1/2Z2 and defining r = x x x 2 , this can be written as 



a = I C( T )&(Ts)dT (9) 

where 



Here 



£{t)= [ ^ r fl-{x 1 )f w+ {^ r )dx 1 (10) 

f w .{x x ) = f 1 -&-{*iMfi(vi)dVi (11) 
J xi yi 

with a similar expression for f^ v+ {x2)- Since the structure functions have been previously 
calculated (see Ref. 8), f^- and can be determined initially, and then only the single 
integral for a above needs to be done. It is important to note that structure functions 
for longitudinal W's differ from those of transverse W's, and thus the scattering cross 
sections must be determined for every combination 
The cross section is given by 

- (12) 
16tts J i _ 4m 2 w /s 

and thus we only need to determine the matrix elements squared for each process. We 
will designate the transverse W's as Wl and Wr, and the longitudinal W as W . 

For the scalar exchange diagram, the only nonzero contributions are WlWl, WrWr 
and Wo Wo- For WlWl and WrWr, the square of the matrix element is 



and for Wo Wo it is 



[s - M 2 S ) 2 + T|M S 2 v * ' \2m 



2 

w 



For the fermion exchange diagram, there are many more spin combinations. The 
full expressions are given in the Appendix. For the vector exchange diagram, the only 



6 



non- vanishing contributions are WlWl (which is identical to WrWr), W£W (which is 
identical to W^Wq) and W Wo. For W L W L , we find 

|yu| 2 , 9_M (l + * 2 )(s-M?)(s-4 m * w )(M? + 2s) 

For ly/VFo", we have the same expression multiplied by s/2m^, and finally, for WoWoj 
we have the same expression multiplied by (s/2m^ + l) 2 . 



5 Results and Discussion 

The results are given in Fig. 7. Results are given for the LHC (y/s = 14 TeV). We see 
that the cross sections are all quite small, never exceeding 10 -3 picobarns. Note that 
we have assumed that the scalar boson flavor changing coupling to the top and charm is 
the geometric mean of the top and charm Yukawa couplings; if it were the top Yukawa 
coupling, the rate could be 0.1 picobarns. Note also that we have scaled the cross-section 
in the case of vector exchange by = 10 -4 , corresponding to the model-dependent 
bound on mixing with the Z; should mixing with the Z be suppressed, this coupling 
could be larger. 

Could such a cross section lead to measureable single-top/single-charm production? 
After all, a 10 -3 picobarn cross section, with an integrated luminosity of 1000 inverse 
femtobarns (possible at the LHC after several years of running), will lead to 1000 events. 
However, there will be a background from single-top/single-bottom production. The 
cross section for pp —>■ qq — > tb, through s-channel W-exchange, is 10 picobarns at the 
LHC|| (the cross section for single top production through gluon-W fusion is larger, but 
the b is much softer and appropriate kinematic cuts can suppress this signal). Thus, one 
would have to be able to distinguish a c in a sample of 10 4 6's. Alhough kinematic 
cuts and good vertex detection would help, and the expected single-top/single-bottom 
cross section can be calculated to a few percent accuracy, it is hard to imagine that an 
unambiguous signal could be detected. It may not be hopeless in the scalar exchange 
case if the coupling of the scalar was unexpectedly large; the rate could then be as high 
as 0.1 picobarns, and one would get a c for every hundred 6's. In this case, the signal 
possibly could be extracted. 

In view of our statements in the introduction, the fact that W + W~ — > ct + tc 
appears to be too small (for the most natural values of couplings) to be detected, while 
t — ► c W + W~ can, in some cases, be detected, may be surprising. After all, the latter 
does have a significant phase space suppression, while the former does not. However, the 
latter also has the possibility of a resonance for some mass range, while this resonance is 
absent in the hadronic collisions. 

We conclude that the process W + W~ — > ct + tc is probably undetectable at the 
LHC. The only hope would occur if the exchanged scalar boson had a surprisingly large 
coupling. The process t — > c W + W~ may be detectable in two cases: (a) With scalar 
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exchange and the scalar mass between 2mw and approximately 200 GeV and (b) with a 
heavy fermion with a mass very near the current lower bound and with an unexpectedly 
high mixing with the second and third generations. 



This research was supported in part by the U.S. DOE contract DC-AC05-84ER40150 
(Jefferson Lab). 



6 Appendix 



Here, we give the results for the fermion exchange diagram in terms of the scattering 

* I 

/I 



angle. Integrating over this angle is straightforward. We define r = cos 2 § and r = sin 2 



2 " ^ — " 2 ' 

and also define K Q = g A \ V t fV*j\ 2 / {t—m 2 F ) 2 . Results are given for all possible combinations 
of a and a' equal to ±1; other cases may be derived from those. 

For WlWl scattering, the square of the matrix element is K s~ 2 (s — M 2 )A, where for 
(a, a') = (1, 1), A = (s - M 2 ) 2 rf 2 (M 2 r + sf); for (a, a') = (1, -1), A = m 2 F f(M 2 r + sf); 
for (a, a') = (—1,1), A is the same with r <-> f, and for (a, a') = (—1,-1), A = 
f(M 2 r + sf)(M 2 f + sr) 2 . For WrWr scattering, the result is the same 

For WrW£ scattering, only the a = a' contributions are nonzero. For (a, a') = (1, 1), 
the square of the matrix element is K Q s~ 2 (s — M 2 ) 3 r 2 f(fM 2 + rs) and for (a, a') = 
(— 1, —1) it is the same with r <-> f. For W^Wr, the result is also the same. 

For WqW f , the square of the matrix element is \K a s~ 2 (s — M 2 )m^ (sf + M 2 r)A, 
where for (a, a') = (1,1), A = s(s - M 2 ) 2 fr 2 ; for (a, a') = (-1,+1), A = s 2 rM 2 F and 
for (a, a') = (—1, —1), A = f(sf + M 2 r) 2 . For the remaining case, (a, a') = (1, —1), the 
result vanishes. Again, for WqW^, the result is the same. 

For W^Wjjr, the square of the matrix element is ^K Q s~ 2 (s — M 2 )m^(sr + M 2 f)A, 
where for (a, a') = (1,1), A = s(s - M 2 ) 2 r 3 ; for (a, a') = (1,-1), A = s 2 fM F and 
for (a, a') = (—1,-1), A = f 2 r. For the remaining case, (a, a') = (—1,1), the result 
vanishes. Again, for WqW f , the result is the same. 

Finally, for M/oW^o, the square of the matrix element is |i^ m^ (s-M 2 )(sf + M 2 r)A. 
For a = a', A = (s — M 2 ) 2 rf 2 and for a = —a', A = sm? F f. 
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FIGURE CAPTIONS 



Fig. 1 This process gives £ — > c W + W~ in the standard model, where F is any 
Q = —1/3 quark. In generic non-standard models, F is a heavy quark. We define the 
upper vertex to be ig^^l + a^ 5 )/2\/2 and the lower to be ig<fi v {\ + ck'ts) /2-\/2. 
Fig. 2 Scalar exchange contribution to £ — > c W + W~ in a multi-scalar model. The 
Yukawa coupling is defined to be —i(a+b'y 5 )g^m c m t /\^2mw and the scalar- vector- vector 
coupling is defined to be igmy/9^vC w . 

Fig. 3 Vector exchange contribution to £ — > c W + W~. The Yukawa coupling is defined 
to be 2g 2 7p(l + «75)/2 and the triple-vector coupling is defined to be gi(g^(k — P) P + 
gup(p — q)[i + gpLi(q — k) u ), where all momenta are defined to flow outward from the vertex. 
Fig. 4 Branching ratio for £ — > c W + W~ due to non-standard model scalar exchange, 
using the vertices from Figure 2. The rate is proportional to C^(a 2 + b 2 ), which we have 
set equal to unity in the figure. 

Fig. 5 Branching ratio for £ — + c W + W~ due to heavy fermion exchange. The solid 
(dashed, dotted) curve is the contribution from 7\ (T 2 ,T 3 ), defined in Eq. 3. Here, we 
have set g\ and g 2 equal to the electroweak gauge coupling, and set \V c f\ and \Vf t \ equal 
to unity; the branching ratio scales with these quantities as shown in Eq. 3. 
Fig. 6 Branching ratio for £ — > c W + W~ due to non-standard model vector exchange, 
using the vertices from Figure 3. In plotting the above, we have set g\g\{\ + a 2 ) equal 
to unity. As discussed in the text, in most models limits from mixing with the standard 
model Z force g\ to be roughly two orders of magnitude smaller than the gauge coupling; 
the branching ratio in such a case will be more than four orders of magnitude smaller 
than shown in the figure. 

Fig. 7 Cross section for W + W~ -> c£ + £c at the LHC (y/s = 14 TeV). For vector 
exchange indicated by the dashed line, the value of has been taken to be 1CT 4 , which 
is the maximum value expected from limits on Z mixing. The solid line corresponds to 
scalar exchange, The other lines indicate the cross section for fermion exchange. The 
labels {LL, LR, RL, RR} refere to the cases where (a, a') are {(—1,-1), (— 1,+1), 
(+1,-1), (+1,+1)} respectively. 
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BR(t --> cWW) 




BR(t --> cWW) 




BR(t --> cWW) 
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